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Abstract 

System and synthetic biology are rapidly evolving systems, but both 
lack tools such as those used in engineering environments to shift the their 
focus from the design of parts (details) to the design of systems (behav¬ 
iors); to aggravate, there are insufficient theoretical justihcations on the 
computational limits of biological systems. To diminish these deficiencies, 
we present theoretical results over the Turing-equivalence of metabolic 
systems, defines rules for translations of algorithms into metabolic P sys¬ 
tems and presents a software tool to assist the task in an automatic way. 


1 Introduction 


There is an increasing interesting, in the academic community, in the modeling 
of existing biological cells as well as synthesis of new ones. Systems biology and 
synthetic biology are parallelly developing themselves with great speed but, in 
some sense, apart and lacking mutual synergy. Nonetheless, a duality between 
these fields may be established and both may benefit of the same set of mathe¬ 
matical, computer science and engineering techniques to improve their analysis 
or design tasks; for instance, molecular computation 16 Section 20.7.5], in 


which the present work is certainly part of, is one of those tools. 

Through the usage of the Church-Turing thesis Section 5.1] p. 183] 
Section 2.5] and algorithmic construction, it will be shown that a particular kind 
of Metabolic P (MP) systems Chapter 3] called positively controlled MP sys¬ 
tems (MPPC) has the same computational power of a register machine and, 
hence, of a Turing machine, along with a demonstration of the equivalence rela¬ 
tion between these models. Then, an outline of the software tool that converts 
a register machine to MPPC is presented. At last, there is a discussion on the 
importance of the results of this work as a connector between systems and syn¬ 
thetic biology and introduction of a computational perspective over these fields 
is stressed. 

For this purpose, the present text is divided as follow: Section and Sec¬ 
tion 1^ introduces the bases of universal computing devices and MP systems; 
Sections and develop the equivalence of the both systems. Section de¬ 
scribes the developed tool and presents an example of translation of register 
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machine to MPPC grammar. Finally, at the last section, we discuss the advan¬ 
tages of the present formalism to both systems and synthetic biology. 


2 Universal Computing Devices 

The concept of (universal) computability is closely tied to concept of algo¬ 
rithm l8|p. 246] that, in turn, is tied with the concept of Turing machine 
p. 246] [19[ Definition 5.17], which may be considered the main computational 
device and stands out by its formal and construction simplicity. 

Although simple, there are other computational devices equivalent to Turing 
machine that are more convenient to use at certain modeling occasions. Some 
examples of equivalent devices are recursive functions, grammars and register 
machine. In the present work, a particular kind of the latter will be used due 
to its simple instruction set and similarity to electronic computer architecture. 


2.1 Register Machines 


According to Minsky 12 Section 11.1], a register machine (or program machine) 
is composed of finite number of registers which has infinity capacity, a small set 
of operations over the registers and a hard coded program, which seems an 
indexed sequence of instructions. At last, an instruction is defined as a triple 
[operation, register reference, list of instructions). 

The basic operations of a register machine must reflect those ones that defines 
recursive function, the (free) movement of the head of the Turing machine across 
the tape, a signal of end-of-computation and limit its operations to the set 
N. Hence, the base register machine model in Table 11.1-1] define four 
operations: zero, successor, decrements or jump and halt. 

[TtI Section 4, p. 225] 


The present work, nevertheless, uses a variation 
p. 293] of the register machine defined by Minsky 12 
seen. 


lill 

Section 14.1] as following 


Definition 1 (Register Machine). A register machine TZ is a computational 
device defined as 

n = [R,0,P) 

where: 


1. R = {Ri, R 2 ,..., Rn} is a finite set of infinite capacity registers, with 
n e N; 

2. O — {INC, DEC, JNZ,HALT} is the set of operations; 

3. P = [Ii, I 2 ,..., In) is the program, with n G N. 

The execution of the program P always start at the first instruction Ii and 
procedures sequentially (unless for programmed execution re-route). 

The instructions I of the register machine TZ are special notations that conve¬ 
niently name operations over addressed registers of TZ according to Definition]^ 

Definition 2 (Instructions). Let the content of register Ri be equal to x. Then, 
it is possible to define the instructions I of the register machine TZ as following: 
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1. INC(Ri) = x + l; 


2. DEC(Ri) 


a; — 1 , if X > 0 

0 , otherwise ’ 


3. JNZ(Ri,Ij) change the execution flow ofTZ setting Ij as the next instruc¬ 
tion to be executed in case x > 0; otherwise, the execution flow keeps 
seguential; 


4-. HALT ends the computation ofTZ. 


A register machine, as seen above, is a simple but powerful computational 
device easily translate to digital computer architecture. Now it is time to ex¬ 
amine that one inspired by cell metabolism, the MP systems. 


3 Metabolic P Systems 


The cell may be seen a small dimension but complex and dynamical system 
limited by a membrane which separates the external world from the internal 
cellular machinery. Acting as an interface, this membrane selectively collects 
molecules from surroundings of the cell and expels others that were produced 
or refused inside it. 

The process of material exchange interfaced by the membrane has a par¬ 
allel with systems’ theory, in which the cell works as a box (isolated system) 
for transformation of substances, suggesting the existence of a computational 
process inside it. Gheorge Paun, aware of this mechanism, proposed a compu¬ 
tational model called P system [14| (precursor of membrane computing) which 
is based on membrane interaction in cell systems but, at the same time, math¬ 
ematically formal and consistent, using the concepts of multiset and rewriting 
systems for the construction of its formal framework. 

Although full of new and interesting ideas for biological modeling, P system 
still presents mechanism too tied to formal languages that forbids its usage to 

p.64]. Attentive 


model real-world metabolic and intra-cellular interaction 10 


to the necessities of bio-modeling, a new membrane computing computational 
model based on this system which is named Metabolic P system or, for short, 
MP system Chapter 3]. 

The primary goal of the MP system is to deterministically model metabolic 
processes, serving as a powerful (discrete) mathematical tool for expressing and 
supporting biological studies in the cell magnifying level; also, it meant to be 
a computational “intermediate language” for easy simulation of the formalized 
models; at last, it should use promptly understandable notation for potential 
users unfamiliar with the theoretical computer science jargon commonly found 
in new computational models. 

Strongly influenced by chemical reactions, MP system has a reaction-like 
notation—that can be seen, also, as a formal grammar one—supported by re¬ 
currence equations and shifts the focus from pointwise string rewriting to a 
population transformation through the usage of conventional mole concept (as 
in chemistry). By the other side, its dynamics is supported by linear algebra 
and relies on matrix operations for solving the recurrence equation system that 
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characterizes the MP system as a (discrete) dynamical one. In technical jargon, 
an MP system can be mathematically represented using the support of a kind of 
formal grammar (named MP grammar). As a grammar object, the MP gram¬ 
mar defines all the rules of an MP system, including the process of multisets 
(through rules and functions), the elements allowed in the multisets and the 
initial state of the systems p. 108]. 

Definition 3 (MP grammar). An MP grammar G is a generative grammar for 
time series defined as 

G = 

where: 

1. M = {xi, X 2 , ■ ■ ■, Xn} the finite set of substances (or metabolites), and 
n € N the quantities of substances. 

2. R = {aj —t /3j I 1 < j < m} the set of rules (or reactions), with aj and 
Pj multisets over M, and m G N the number of reactions. 

3. I = (a;i[0],a;2[0],... ,a:„[0]) is the vector of initial values of substances or 
the metabolic state at initial step (step 0). 

4 . ^ = {(pi, p 2 , ■ ■ ■, Pm} is a set of functions (also called regulators), in which 
every pj : R" 1 —>■ M, for I < j < m, is associated with a rule rj G R. 

According to definition]^ G generates (a set of) time series, each of them 
representing the “amount of quantity” of the substances during the time and its 

time series is calculated for any time t if and only if — € N for a given constant 

T 

T. 

Notwithstanding, the rules Ui ^ Pi G R depends, as equivalently happens in 
chemical reactions, on the quantities of the “substances” in the system, which 
can be expressed with the support of two different concepts: one that maps 
the multiplicity expressed in the rule for a substance to the actual number of 
molecules (of the substance) in the system, and; the quantity of mass the unit 
of the multiplicity represents (for a particular substance). 

Hence, if the aforementioned three restrictions are that in consideration 
along with an MP grammar G, it is formally defines a MP system. 

Definition 4 (MP system). A MP system M is a discrete dynamical defined 
as 

M = (G, r, V, p) 

with 

1. G being an MP grammar following the definition^- 

2. T S R, the period (amount of time) of a computational step; 

3. 1 / G R, the number of molecules that represents the (conventional) mole in 
the system; 

4 . p € R" is the vector of the mole masses of substances. 
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As described in definition the quantities of the substances are dependent 
of its previous values and a variational function that may depend on other 
substances and parameters (in the case of the parametric MP system). This 
additional variance through time is represented as a recurrent equation which 
the future value of a substance X is represented as a;[i + 1] cx x[i\. 

For the computation of these step values, nonetheless, two mathematical 
accessories were developed. The hrst, the stoichiometric matrix, is based on the 
arithmetic executed over chemical reactions to calculate the balance of molecules 
in a chemical system; the other, equational metabolic algorithm, synthesizes the 
whole computational process specified by the an MP system. 

Definition 5 (Stoichiometric matrix). Let = Oi —>■ j3i, where on (with an 
equivalent for /3i) is represented as X Xj I kij GN AXj G M. 

Let mult^(Xj, Tj) = fcA be the multiplicity, for the right side (at) of the rule 
Vi, of the substance Xj in the rule. Similarly, there is mult~(Alj, r^) = k~j for 
the left side (Pi) of the rule. 

A stoichiometric matrix A, of dimension \M\ x |i?|, has each of its elements 
defined by 

apra = mulf'‘(X/,rm) - mult“(X/,rm) 
with ,1 < I < \M\ and 1 < m < |i?|. 

Definition 6 (Equational metabolic algorithm—EMA). Let U[i] = {(pi (i), <p 2 {i), ■ ■ ■, Pm (*))^ 
be the vector of values, in the time step i, of all regulators, and A the stoichio¬ 
metric matrix. 

The vector of substance variation at step i, A[i], is computed by the equation 

A[f] = A X U\i] 

so-called Equational Metabolic Algorithm whom computes the value of any sub¬ 
stance in the time future time step i + 1 through the recurrent equation 

X[i + 1] = X[i] + A[f] 


The above definitions, now, complete specihes the discrete dynamical system 
and the ways to compute its values. 

3.1 Positively Controlled Metabolic P System 

The cell, motto of the MP systems, keep itself working based on metabolic rules. 
Those, of biochemical equations nature, require enough quantity the metabolites 
on the left-hand side of their equations to transform the matter and, hence, keep 
the metabolic circuitry active. Thus, if enough quantity of all metabolites are 
provided for all equations, the metabolic circuitry works properly; in the total 
absence, it does not work. However, in the case which a subset of rules do not 
have enough quantity of some metabolites, these rules are inactivated while the 
others keep their activities normally. As a result, the metabolic circuitry of the 
cell keeps working with a broken chain of reactions. 

This kind of operation on cell metabolism happens because the (mathemat¬ 
ical) operations in cells (and bio-chemistry, in general) happens in the space of 
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natural numbers, which implies that half (decimal) quantities or debt (negative) 
quantities cannot be considered in the (left-hand side of the) rules. 

In the MP systems world, the above cell-like behavior is defined in a subclass 
called Positively Controlled Metabolic P systems, or MPPC: it is a standard 
MP systems with few controls to ensure all operations are executed in the set 
of positive numbers. 


Definition 7 (MPPC Grammar). A MPPC grammar Q+ is a standard MP 
grammar Q respecting the following restrictions, at each computational step Si: 


P\s, 



if K> 0 
otherwise’ 


for all (fi G d); 


2. ^ < X, with $ 2 , the set of all fluxes related to rules in which the 

metabolite x is in the left-hand side of the rule; otherwise, (p = 0, Wip G d>“ 
at the execution step. 


As it will be seen later in the next sections, this definition is useful not 
only for cell modeling, but also for establishing a comparison with universally 
computational devices. 


4 Computational Architecture as Metabolic Sys¬ 
tems 

One of the goals of synthetic biology is to produce programmable metabolic 
systems that could, just as a designed machine, perform well-defined tasks for 
a certain objective. Some successful approaches of computer-aided design of 
biological systems already exists [^, but they are limited by the expressiveness 
of digital gates. 

Theoretically, nonetheless, we are possible to synthesize metabolic systems 
more significant (in computational power) than those of Boolean circuits. For 
the present purpose, a programmable and Turing-potent (register) machine is 
implemented in a metabolic P system. 

Theorem 1 (Translation of Register Machine to Positively Controlled MP 
grammar). For any register machine TZ exists an equivalent positively controlled 
MP grammar . 

Proof of Theorem^ Given a register machine TZ = {R,I,P) with |i?| = r and 
|P| = p, a positively controlled MP grammar Q+ = (M, Ru, /, <&) is constructed 

1. adding a metabolite Ri in the set M for each register Ri G R; 

2. adding a metabolite Ij in the set M for each of the instructions in Ij G P; 

3. adding a metabolite Lj in the set M for each instruction Ij G P oi the 
type JNZ; 

4. adding a HALT metabolite in the set M; 
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5. defining the initial state of the metabolites Rj equal to the initial values 
of the registers Rj, the initial values of all the other metabolites to 0 and 
the initial value of Ii to 1; 

6. adding the rules to Ru and the fluxes to $ according to the following rules: 

(a) if Ij is INC or DEC, then Ij Ij+i : Ij] 

(b) if Ij is INC(Ri), then % ^ Ri : Ij] 

(c) if Ij is DEC(Ri), then Ri ^ : Ij] 

(d) if Ij is HALT, then Ij —>■ HALT : Ij] 

(e) if Ij is JNZ(Ri, lie), then 

i. Ij —>■ Lj : Ij] 

ii. Lj —>■ Ik '■ Lj — Ij+i] 

iii. Lj ^ Ij+i] and, 

iv. 0 — Ij-\-i ■ Ij — Ri- 

From the rules above, it is possible to notice that Ij and Lj instructions 
controls the execution flow of the system and satisfies 

HALT + = 1 






L,< 


ensuring no two instructions are executed at the same time, but its execution 
starts from instruction Ii and proceeds sequentially (or jumps to another one 
in case of a satisfying JNZ instruction). 

All operations are mappings from and to the N set once both R. and , by 
definition, restrict their operations to this set. 

At last, when a rule Ij —)■ HALT is performed, the system is stuck in a fixed 
point since there is no rules for “exiting” this state. □ 


5 Translation Strategy 

Theorem [l] defines the relation between register machine and positively con¬ 
trolled metabolic P system, but it does not develop the intuition for its reasons. 
For a complete comprehension of this equivalence and, consequently, its im¬ 
pact in correlated fields, this section will focus on the constructive analysis of 
a MPPC from a register machine (just as proposed by Theorem and discuss 
the challenges of this procedure. 

5.1 Is It Possible To Translate Register Machines Into MP 
Systems? 

A register machine is a general purpose computational device presenting the 
same computational power of a Turing machine; it allows the implementation 
of any algorithm and can deal with languages in the whole spectrum of the 
Chomsky hierarchy p. 272]. By the other side, MP systems have been show 
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to be equivalent to Petri networks and have been used to model complex sys¬ 
tems, from (originally) metabolic systems to mathematical series Chapter 3]; 
however, MP systems are not proved to be as powerful as a register machine 
and an equivalence between these systems cannot be set yet. 

Nevertheless, MP systems presents a series of characteristics that induces to 
the idea of an existing equivalence between both systems. It is important to 
note that these features, by themselves, do not ensure the existence of the any 
kind of relation between the formalisms, but solely encourages the investigation 
through a number of theoretical arguments. 

At first, MP systems are (a particular kind of feedback) dynamical systems] 
therefore, it receives a set of states, manipulates them and feeds them back to 
itself in order to keep the computation going on. It is a computational device 
with proved restricted power [^, but with potential for increasing it 0 m 
particularly because it does not restrict the set of functions it may works on— 
i.e., the fluxes may be any computable function according to its definition. 

Then, both systems works in a discrete representation of time, computing 
the elements in well-defined steps instead of continuously. Each step executes a 
finite number of instructions in a finite amount of time and, more, with a finite 
representation (i.e., there is no guarantee of infinite precision for every object 
computed). Although these limitations on time are not strong enough to ensure 
the equivalence (after all, both Turing machines and finite state machines are 
discrete but not equivalent), it prevents analyzes that goes beyond the Turing 
equivalence, named super-Turing ID (I- 

At last, MP systems have representations with origins on (or inspired by) 
formal languages (grammar), rewriting systems and graph theory (MP graph 
p. 109 and Figure 3.1]). Not enough, it is part of membrane computing, a 
wide field with other devices as powerful as (or even more than) the Turing 
machine p. 179] (with properties, though, that make them impracticable for 
daily basis use). 

5.2 What Is The Strategy To Implement? 

In order to show that MP systems are general purpose computational devices 
as powerful as register machines (and, hence, Turing machines), it is enough 
to show that it is possible to implement in MP systems the same algorithms 
accepted by the register machine. And to satisfy the generality of these imple¬ 
mentations, it suffices to translate the register machine instructions into MP 
equivalents and show they work properly for a well-defined subset of all accept¬ 
able algorithms in that machine. 

Expressions of behavior in MP systems are described solely through its MP 
grammar Definition 3.1], which is defined as a quadruple consisting of sets of 
metabolites, rules, initial states and regulators (or fluxes). Therefore, the defi¬ 
nition of the four elements of a MP grammar defines univocally the operational 
actions of a system; if the definition is relaxed and the initial states are not 
specified, the relaxed MP grammar defines a whole class of behaviors instead of 
a specific one. 

The strategy, then, is to restrict the definition of a relaxed MP grammar that 
models the whole class of problems accepted by a register machine^ populating 
the sets that compose it according to the next strategies. 


5.2.1 Metabolites Set 


The metabolites set M in a MP system is the equivalent to M = X U F, the 
union of the state space X and output value space Y of dynamical systems |^. 
It declares all the variables (states) of the system, those allowed to change its 
values through the computation process. 

Analyzing the register machine instructions 17 Section 4, p. 225] Chap¬ 
ter 11] from Section (the unique procedures that alters the states of a register 
machine), it is possible to recognize and isolate all the variables of the machine 
just looking to the addressed properties subject to manipulation. These can be 
classified as 


• Registers. Three out of the four instructions (INC, DEC and JNZ) refers 
to registers, either for manipulation or query the data. Registers are the 
external variables pp. 74] of the register machine, the ones of interest 
for an user that feeds the machine with a program. 

• Instruction Pointer. Each of the instructions of the machine are executed 
sequentially, one at each computational step; however, this serial behavior 
can be modified with the usage of the JNZ instruction, which “jumps” 
the next execution step to another indexed instruction if an addressed 
register has its contents different from zero. And in order to keep the 
track of the instructions to be performed, there is a special kind of register 
(called program counter) whose content is the index of the instruction to 
be performed. Normally, its initial value is one and it is automatically 
incremented by one at the end of the execution of the instruction; however, 
JNZ is able to change it to whatever (valid) index it is necessary, hence 
qualifying program counter as a register machine’s variable. 

• Final State. To indicate the end of a computational process executed 
by a register machine, there is a special instruction called HALT which 
signals the machine to stop its execution because the algorithm in process 
has finished its operation. Hence, halt {final or acceptance state [^) is a 
particular that marks the end of the computational process. 


5.2.2 Rules Set 

The rules in a MP systems define the relation between the diverse variables of 
the system, without pay attention to the value update of the variables. Rules, 
then, characterizes the operation of the system and are presented in a —>■ /3 
structure—which can be read as a becomes (or transforms itself into) ft. 

In the current scope, the definition of the structural behavior of the general 
MP system that simulates a register machine is done through the observation of 
both the register machine and the metabolite set: the former is detailed studied 
through the standpoint of the instruction set and the possible program flows 
(flowchart) and describes the behavior being reproduced in the new system; 
the latter provides the description of all the systems’ states (variables), the 
principal constituent of the rules and defines the three subdivisions of behavior 
of the register machine: manipulation of the registers, manipulation of the flux 
of the instructions ’ execution and final state of the computational process. 
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Manipulation of the Registers The INC and DEC instructions are dedicated 
to manipulation of the registers, while JNZ uniquely consults their values. 

The instructions that alters the values of the registers, consequently, are 
those to induce changes in (some of) the system’s states (the external ones 
p. 74]), expressing one of the (main) behavior of the system; therefore, they are 
modeled in the MP system through particular but very simple rules: they just 
add to or remove values from addressed registers. 


Table 1: MP rules for register machine instructions that manipulate registers’ 
values. _ 


Register Machine Instruction 

MP Rule 

INC(Ri) 

0 ^ 

DEC(Ri) 

Ri^% 


The first of the rules of Table can be interpreted as “from somewhere the 
register i?* receives a value” (and the second one as “the value of the register 
goes to somewhere”) because of the structure of MP rules: inside a delimited 
environment {e.g., a cell), the principle of mass conservation must be preserved 
and any addition (or removal) of metabolites quantities that breaks this principle 
must “come from (or goes to) somewhere”; somewhere, in this context, must be 
understood as outside the membrane delimiting the MP system. 

Manipulation of the Flux of Instructions Not all the states of system 
are registers; those responsible for the sequential execution of the system and 
the control of its computational workflow, the instruction pointers, represents 
the computational machinery and hide the trickiest translation process between 
register machine instructions and MP rules. The evaluation of these variables 
may be divided in two different subclasses: sequential execution of instructions 
and manipulation of the execution’s workflow. 

The first of the subclasses, sequential execution of instructions, consists of 
situations when the token of active instruction is passed from one instruction 
to another consecutive. Thus, let / = {Ji, / 2 , • • ■} be the set of indexed instruc¬ 
tions of a register machine representing an algorithm. The computation of the 
algorithm represented by / is defined as the sequential execution Ji h /2 h ..., if 
and only if Ij is a instruction of the type INC or DEcQThen, from the previous 
statement, it is trivial to define the following conversion rule: 

Let be the and (j -I- 1)*^ indexed instructions, with 

Ij belonging to one of the types INC or DEC. Then, there is a 
MP rule Ij Ij+i that sequentially executes the instruction Ij 
in a MP system. 

For the above rule, it is important to notice that there is the restriction of 
the type of instruction (INC or DEC) solely for the Ij instruction; the Ij+i, on 
the other hand, can be any of the register machine instructions. 

^Although JNZ instruction also pass the token in a sequential manner, it is not true for 
every performance of it; therefore, this instruction is excluded of this subclass and is analyzed 
later. 
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The other subclass, manipulation of the execution’s workflow, comprises the 
situations in which the computation of / does not follow the sequential disposi¬ 
tion of the instructions, but changes the execution order under certain circum¬ 
stances. For this, it is required the JNZ instruction, the exclusive instruction 
capable of performing such re-arrangements. 

Dissecting the JNZ instruction (Definition]^, it is possible to see it has two 
arguments, a register address and an instruction address, one used to perform an 
inequality comparison against zero of a register value and the other to redirect 
the current workflow in case the comparison results are true. be.if the addressed 
register Ri does not have its value equal to zero; since this instruction is inher¬ 
ently composed of two operations, it persuades us to convert it in, at least, two 
MP rules. 

Figure 1: Flowchart of the JNZ(Ri,Iij) instruction. 



At first, let us focus in the comparison operation of JNZ. Looking at Figure]^ 
it is easy to see the detour of the sequential instructions’ execution happens when 
the value of the register Ri differs from zero; this remark suggests most of the 
trickery will be developed in order to make this execution branch feasible in MP 
systems. 

The artifice is to create a metabolite Lj that represents the virtual com¬ 
parison instruction of instruction Ij = JNZ(Ri, Iji) and, then, let it support the 
control of the redirection of the instructions’ execution. 

Thus, Lj must receive the token previously in Ij and, in your own computa¬ 
tional step, choose whether Ij+i or R will be the next instruction to be carried 
out. 

For the reception of the token, it follows the same mechanism previously 
used for sequential execution of the instructions: the Ij token becomes the Lj 
token with the simple rule Ij —>■ Lj. 

The choice of the instruction path, in turn, makes extensive use of regulators 
(which will be discussed ahead in Section 5.2.3) in order to produce its rules. 
Intuitively, though, it is possible to develop them. 

When the contents of 7 ^ 0, Lj must redirect the execution flow to the 
instruction R through the transfer of its token to this new instruction; therefore, 
again, the sequential rule style operates and the produced rule is Lj —>■ R. 

However, when = 0, a problem arises: neither R —>■ R+i nor Lj —>■ 
Ijj-i rules may be used to give, in the context of a metabolite becomes another 
metabolite, Ij+i the execution token—the former {Ij) does not have tokens 
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anymore and the latter (that would be the correct rule) was used to transform 
the token from Lj to Ik- The solution, then, is to create two rules to reproduce 
the behavior of Lj — >■ Ij+i'- one takes the token from Lj and the other gives it 
to Ij+l- 


L^ —y I- 


j+i 


Lj 


ij + 1 


At end, an instruction Ij = JNZ(Ri, Ijt) becomes four MP rules summarized 
in Table m 


Table 2: MP rules for the Ij = 

JNZ(Ri,Ik). 

Functional Description 

MP Rule 

Pass token to pointer of virtual 

^3 

comparison instruction Lj 

If Ri 7 ^ 0, pass token to Ik 

Lj Ik 

If Ri = 0, pass token to 

Lj —>■ 0 

0 —)■ Ij+i 


Final State of the Computational Process At the end of a computational 
process, some kind of flag is necessary to signalize its termination, otherwise 
there is no way to guarantee we are seeing a partial result of the computation 
or the system’s hnal state. Turing machines, in their formal dehnition Defi¬ 
nition 4.1.1], make this explicit signalization through the halting states set, and 
the register machine uses its HALT instruction; these halting objects, as the name 
suggests, stops the computational device execution. 

By the other hand, MP systems are feedback discrete dynamical systems 
and the concept of stopping the computation does not exist in this context. In 
other to imitate this behavior, it is possible to add to the systems’ variables 
an additional state (which we will conveniently call it halt state) which marks 
the end of the systems’ computational process and, at the same time, is a fixed 
point of the dynamical system. For this purpose, it is enough the addition of a 
single halt state that is the destination of the token from the HALT instructions 
and the absence of rules with the halt state in the left-hand side of it. Hence, 

Let Ij be a HALT instruction. Then, there is a MP rule Ij —)■ 

HALT that models the end of the computational process with 
signalization of the final state. 

Clarifying the above definition, a register machine does not restricted your 
programs to have an unique HALT instruction and, thus, it may have several of 
them, all of them going to a single halt state. Therefore, the “token” of the j**' 
instruction becomes the signal of end of computation and no other instruction 
is executed, because: (i) the value (“token”) that controls the instruction flow 
is passed from the instructions (Ij) to the halt state (HALT) and (ii) there is 
no rules in which the halt state is consumed {i.e., in the left-hand side). 
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In fact, those two arguments also justify the halt of the device as a fixed 
point of the dynamical systems: with no rules being executed, no modification 
in the states are performed and, that being so, any other execution step of the 
MP system produces the same, previous state {Shalt b Shalt b ... b Shalt)- 

5.2.3 Regulators Set 

In a MP system, the metabolites set defines the state variables of the system. 
The rules set, its skeleton behavior. But the update of the values of the state 
variables depends on mathematical functions associated with the rules, also 
known as regulators. 

The regulators computes the rates in which the metabolites in the left-hand 
side of rules are transformed into the ones in the right-hand side. In both 
mathematics and biology, regulators are often called flux. 

Since does not make sense to study regulators detached from rules, the clas¬ 
sification of the equivalent regulators of a register machine will follow similar 
subdivision used for the rules: update of the registers’ values, update of the 
instruction pointers and signalization of end of computation. 

Update of the Registers’ Values The rules presented in the Table de¬ 
scribes the process of updating a register value, but not the functions which 
actualizes it; these functions (the regulators) can be easily inferred from the 
specification of the register machine instructions. 

Initially, the INC operation: it is a primitive recursive function which adds a 
single unit to the actual value, or INC(Ri) = Ri + 1, which can also be stated, in 
a regulator perspective, as “INC(Ri) changes the value of Ri with the constant 
rate of 1”. The latter declaration, in turn, denotes the value of regulator as 
the constant 1; then, it is easy to define the regular for the INC instruction: 
fine = 1 - 

Conversely, DEC follows a similar reasoning, DEC(Ri) = — ][^ which can 

also be stated as “INC(Ri) changes the value of Ri with the constant rate of 
— 1”. Nonetheless, it is also necessary to notice the DEC rule “consumes” the Ri 
value is (because Ri is in the left-hand side of Ri 0) and, hence, its flux is 
already negative. Thus, :/Jdec = 1- 

Update of Instruction Pointers and Final State As with the case of 
the update of the registers’ values, the behavior of the instruction pointers are 
well defined in the Section [5. 2.2 [ but no modification to the state representing 
the instruction pointers are yet defined or performed. But before defining the 
regulators for their rules, it is important to reflect a little about the nature of 
the instruction pointers. 

The instruction pointers states are flags that, when activated, indicates 
which of the instructions are being processed in the current computational step; 
they are the stratagem used to transform the parallel execution of rules in MP 
system (matrix multiplication at Definition into sequential ones as present in 
the register machines. 

^Actually, DEC{Ri) = maxfiti — 1,0} since it is an operation N i—> N. It is guaranteed by 
the positive control of MPPC systems and will be discussed, again, further in this section. 
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The idea of the instruction pointers are similar to the program counter of 
counter and random access machines, with the particularity that each instruc¬ 
tion has its own instruction pointer, they are limited to binary values (deacti¬ 
vated and activated, respectively 0 and 1) and no two instruction pointers may 
be activated at the same time, formally defined a^ifaZt -|- = 1. 

Thus, it is easy to define the regulators of the rules for sequential execution 
of instructions: for each rule of the form Ij —> Ij+i or Ij —> Halt, the regulator 
is defined as ippointer = Ij- 

(The Halt state variable may be seen as a particular case of the instruction 
pointer ones: it marks the end of the computational process instead of the 
instruction being computed.) 

Finally, we must carefully design the regulators for the previously defined 
Ij = JNZ(Ri, Ik) rules. Thus, let us recall the mechanism of this function: 

it compares the value of the register Ri against the number zero 
and if it is equal, the next instruction to be executed will be 
the sequential instruction Ij+i, otherwise it will re-route the 
execution flow to the instruction Ik- 

At first, then, JNZ delegates to the virtual comparison instruction Lj the 
comparison of the register Ri value against zero. Modeled by the rule Ij — >■ Lj, 
it is actually a particular kind of sequential instruction execution, but inside 
the JNZ rule. For this reason, its regulator follows the previous defined for the 
sequential execution rules and assumes the value of Ij. 

Next, there are two rules for the case in which Ri = 0: one passes the token 
to the sequential instruction Ij+i and the other certifies the redirection to the 
Ik instruction will not ever occur. The former case, represented by the rule 
0 —>■ Ij+i, is a variation of the sequential execution when Ri = 0, what drives us 
to consider Ij (sequential execution) and Ri (register value) as elements for the 
regulator. Observing the Table of the possible values of these metabolites, it 
is easy to infer its regulator as = max(/j — Ri,0). 


Table 3: Values of Ij+i depending on the values of Ij = JNZ(Ri, 1^) and Ri. 


Ij 

R, 

Ij+1 

0 

0 

0 

0 

7^0 

0 

1 

0 

1 

1 

7^0 

0 


For the latter of the rules when Ri = 0, Lj —>• 0, its regulator just have 
to control that the virtual comparison instruction Lj will not be active (be., 
its value will be zero) when Ij+i is; given that Lj is solely activated by the 
instruction by a sequential-like rule, it is certainty it will not get any other 
increase of value outside the JNZ workflow and, hence, it has a maximum value 
of one. Then, it is trivial to infer that = Ij+i, establishing Ij+i as the 

repressor metabolite of Lj. 

®The presence of Halt is justified by the necessity of the system be in execution mode, not 
halted. 
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Finally, the lasting rule of JNZ concerns to the re-route of the execution flow 
to the fc**' instruction of the algorithm, lu. As previously state, a redirection of 
the flux of execution happens when the value of the register Ri is different of 
zero, i.e., it is a “secondary level” task inside the JNZ workflow, which depends 
on the rules when Ri = 0. By the other hand, the redirection happens, as state 
by its rule Lj —, as a consequence of the activation of the virtual comparison 
instruction. Minimizing the metabolites set that influences the re-route sub-task 
of JNZ and combining all their possible values in the Table the pattern that 
arises between this table and Table [^correctly induces to conclude the regulator 
for this last rule is — max{Lj — Jj+i,0}. 


Table 4: Values of Ik depending on the values of Lj and Ij+i- 


Lj 

Ij+1 

h 

0 

0 

0 

0 

1 

0 

1 

0 

1 

1 

1 

0 


It is interesting to notice, however, that both 0 Ij+i and Lj R have 
regulators of the form max{X — V, 0}, or subtraction in the natural set of the 
numbers, in order to restrict the regulators to become negative numbers. This 
guard, also present in the DEC instruction, provides two curious features to the 
system; 

1. the comparison operator {greater-than) required in JNZ, ensuring a posi¬ 
tive, not null regulator if and only if A > V; and, 

2. preservation of the coherence of the rules. 

While the former characteristic is evident and does not require further explana¬ 
tion, the latter preserves the structural rules of the system: a negative regulator 
inverses the behavior of the associated rule, i.e. 

{X :-e) = (Y ^ X : e) 

Hence, the concession of negative fluxes in the system would allow the modifi¬ 
cation of its set of the rules, with the modihcation of a certain rule r^ = {X —>■ 
Y,Y —t A} depending of its regulator—or, in other words, the system could 
represent two different programs depending of the positiveness of the regular of 

Ti- 

Furthermore, this restriction pattern goes in accordance with two other (and 
important) principles: the positive control of MP grammars, which asserts 
that fluxes must be greater or equal to zero, and computability over the natural 
numbers Definition 4.4.1, Figure 4-19 and Section 4.7] Section 2]. 


5.2.4 Program and the Initial States 


In Section 5.2 we introduced the concept of relaxed MP grammar as a variation 
of the MP grammar definition without the requirements to set the initial states. 
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This flexibility allowed to construct the equivalence instruction of register ma¬ 
chine in MP as general as it is possible, without restraining the systems to a 
single, particular behavior: a relaxed grammar defines a whole class of dynamics 
based on general rules, not on specific simulations dependent of initial states 
values. 

Although it is enough to define an equivalence with register (and Turing) 
machines through bijection of its composing elements, relaxed MP grammar 
does not allow the definition of the concept of program pp. 210-211 and 
Definition 4.4.1] [12[ p. 202] because nothing such as initial configuration 
Definition 4.4.2] is specified. 

The initial configuration solely defines values for the states of the system and 
no modification of the designed model is undergone. It acts as a constrainer of 
the device model to ensure the correct execution of the dynamics. 

Reviewing the states (metabolites) set of the MP system equivalent to the 
register machine. Section [5.2.1| and later addition in Section [5.2.2[ it is possi¬ 
ble to partition them under two categories: (i) the memory units, which com¬ 
prise the registers-equivalent metabolites (i?i); and, (ii) the instruction pointers, 
containing the instruction pointers Ij, the halting state Halt and the virtual 
comparison instruction Lj. The first kind is, as in similar computational mod¬ 
els Section 2.1.1] Section 11.1] Section 4.4], memory units that keeps 
the state of the internal computations but also serves as an interface for input 
(and, at the end of the computation process, output) data in the program; thus, 
it must always be initialized with zero values when representing internal compu¬ 
tation states, while is open to be freely set to (positive) values when representing 
storage for input data. 

The second kind, instruction pointers, are internal state that control the exe¬ 
cution of the MP system as a sequential, computational device. The instruction 
pointers and virtual comparison instructions, together, represent the program 
counter Defintion 4.4.1] in the MP formalism while the halt state indicates 
the end of the computation. 

In order to correctly execute a program from a register machine, nonetheless, 
it is also necessary to specify the states of the instruction pointers satisfying the 
two following constraints: 


h = 1 ( 1 ) 

|K| 

Ki = 1, for any t gN and Ki G K. = Ij U Lj U Halt (2) 

i=l 


Equation § defines a sequence S = (/i, / 2 ,...) for the instructions codified 
as MP rules merely with the definition of its first element, while the other 
sequential elements follows from the sequential execution rules (Section 5.2.21. 
Equation ([^, in turn, restricts the execution to solely one instruction (pointer) 
at each step, guarantees the device is either on work or halted and, finally, forces 
all the states in the instruction pointer category to have initial values equal to 
zero (except for Ii = 1). 
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6 Software 

Following the theoretical development that derived the Theorem a piece of 
software were codified in order to provide experimental evidence of the validity 
of the theory, as well as a tool for automatic translation of register machines 
specifications into MP systems—which, later, expanded for also simulated both 
systems and ensure their equivalence. Here, though, the bisimulation is not ver¬ 
ified, but instead a verification of the final state of both systems is performed}^ 
However, if one of the systems does not halt, it is not possible to perform the 
comparison of the systems, a expect situation and in accordance to the halting 
problem Section 5.3] Section 5.1]. 

The existing software, a command-line tool available for the major operat¬ 
ing systems players, can be divided in three main parts: translator^ executor 
and comparator, as represented in Figure]^ This modular structure allows an 
unique flexible feature to the tool to choose the atomic operations to perform 
for a determined task; also, it provides the opportunity to implement parallel 
execution of modules and easy addition of new components. 


Figure 2: Dataflow of the software. 

Register 

Machine 

Code 


i 


Translator 


Memoi^ 


Representation 


Memoi^ 


Representation 



Initial ^ 

Executor 


Executor 

State ^ 

Register Machine 


MPPC 



Initial 

state 


When called by the command-line, the software may solely translate the 
register machine specification into a valid MP grammar (Figure]^ or generate 

■^The bisimulation is not necessary to be performed in the software due to the result of the 
Theorem ^ which states an equivalence. There is space, though, for this verification if the 
intermediate steps of both simulations are compared according to [151 Definition 1.4.2]. 
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the equivalent MPPC system to the register machine, simulate both of them and 
print if they are equivalent based in two criteria: both system halted and the 
final state of the registers and their respective metabolites in the MP grammar 
have the same values (Figure [^. 


Figure 3: Register machine specification for the max(i?i,i? 2 ) function. 

JNZ(1, 4) 

INC(4) 

HALT 

JNZ(2, 7) 

INC(3) 

HALT 
INC(5) 

DEC(l) 

DEC(2) 

JNZ(5, 1) 


Figure 4: Translation from register machine to MP grammar with the command¬ 
line topl^_ 


Is . /fegtomp.exe -i esempi/max/max.txt 
IConpilation: [SUCCEED] 

Ho -> R4 : 12 



19 -> HO : 19 

110 -> 0 : no 

0 -> LIO : no 
LIO -> 11 : LIO 


Figure 5: Simulation and comparison of register machine and equivalent MP 
grammar with the command-line tool. 


S ./regtomp.exe -i esempi/max/max.txt --state-counter esempi/max/regs.txt --state-mp esempi/raax/mp.txt 
[SUCCEED] The register machine and the MP system are equivalent. 

Registers: [(1,2),(2,0),(3,1).(4.0).(5,3)] 

MP: [2,0,1,0,3,1,0,0,0,0,0,0,0,0,0,0,0,0,0] 


In the previous example, both devices perform the computation (i?i, i? 2 , Rs, I? 4 , Rs) 
* 

(5, 3,0,0,0) h (2,0,1,0,3) equally, writing 1 at to indicate the value of the 
register i?i > i ?2 (otherwise, i ?4 = 1 to indicate Ri < i? 2 )- 
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7 Conclusion 


Among the several models of cells and metabolism, it is difficult to find one 
that naturally models the cellular system and provides significant insight to 
challenge the biological machinery against existing and well-known synthesis 
theory, such as digital circuits or computer programs. Thus, these models are, 
usually, very useful for the systems biology analysis, but they lack contribution 
for the counterpart synthetic biology. 

The present work extends the existing Metabolic P system, extensively used 
for analyzing biological behavior, to computation theory and provides theoret¬ 
ical and software support to convert programmable instructions into positively 
controlled Metabolic P grammar. Particularly, it stands out from competing 
models because (i) MP grammar is the formalism for every metabolic process; 

(ii) it is a proven Turing-powerful model (under positively controlled restriction); 

(iii) MPPC mimics the cell behavior in the borderline situations; (iv) it provides 
software tools to transform software specification into MPPC description. 

These features, when combined with hardware synthesis techniques 13 or 
database of metabolites (similar to [^), easily brings the ambition of lab-on- 
chip and synthesis biology to the laboratories day-by-day reality. 
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